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Abstract 

Due to its large surface area and strongly attractive potential, a 
bundle of carbon nanotubes is an ideal substrate material for gas stor- 
age. In addition, adsorption in nanotubes can be exploited in order to 
separate the components of a mixture. In this paper, we investigate 
the preferential adsorption of Z?2 versus H2 (isotope selectivity) and 
of ortho versus para (spin selectivity) molecules confined in the one- 
dimensional grooves and interstitial channels of carbon nanotube bun- 
dles. We perform selectivity calculations in the low coverage regime, 
neglecting interactions between adsorbate molecules. We find substan- 
tial spin selectivity for a range of temperatures up to 100 K, and even 
greater isotope selectivity for an extended range of temperatures, up 
to 300 K. This isotope selectivity is consistent with recent experimen- 
tal data, which exhibit a large difference between the isosteric heats 
of D2 and H2 adsorbed in these bundles. 



1 Introduction 

The adsorption of gases within single-walled, carbon nanotubes (SWNT) 
has recently attracted broad attention among physicists, chemists, materials 
scientists and engineers H g, |, § |, g, % |, g 0, Q 13 S S EI 
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Experiments have shown that one can create ordered arrays of nanotubes 
H, |19|, m, H, which form a close-packed bundle (or rope) of SWNT's. 



Ideal bundles of nanotubes consist of very long strands of nearly parallel 
tubes held together at equilibrium separation by intermolecular forces. The 
small diameter and large aspect ratio of the nanotubes make them interesting 
systems for gas adsorption. In such a bundle, physisorption may occur at four 
distinct sites: i) inside the tubes (endohedral adsorption); ii) in the interstitial 
channels (IC) between three contiguous nanotubes (exohedral adsorption); 
hi) in the grooves between adjacent nanotubes on the external surface of 
the bundle, and iv) elsewhere on the outside surface of the bundle. In this 
paper, we assess the ability of nanotube bundles to preferentially adsorb 
specific isotope and spin species. 

Small atoms and molecules, in contrast to larger molecules, have been 
predicted to fit well in the ICs of a nanotube bundle ITR E31. For such 



species, the ICs environment provides a large number of neighboring C atoms 
at nearly optimal distance from the adatom, so the interstitial binding energy 



is larger than is found in other known environments (16|. The grooves also 
present strongly attractive environments for both small and large molecules. 
The binding energies in these sites have been determined experimentally to be 
between 50 % and 100 % greater than the values on graphite j|0[ |42|, [43|, |37j . 

Molecular hydrogen inside SWNTs is a particularly appealing system to 
study, since its adsorption in ICs and grooves is relevant to both gas storage 
and isotope separation |2B|, [ft], [JTJ. "Selectivity" is the term used to 
describe the separation, or selective adsorption, of one species relative to the 
other species of a mixture. In a mixture of two components, the selectivity 
of component 1 relative to that of 2 is: 

S = {x 1 /x 2 )/{y 1 /y 2 ) (1) 

with Xi(yi) the pore (bulk) molar fractions. Since two isotopes have similar 
sizes, shapes and interaction potentials, the separation of isotopic mixtures is 
a difficult and energy-expensive process, requiring special experimental tech- 
niques, such as cryogenic distillation, diffusion separation, laser isotope sepa- 



ration or microwave molecular separation p7|, H) HI- However, most of these 
processes have low selectivity for separating H 2 isotopes. Recently, a novel 
separation technique called quantum sieving was predicted to be particularly 



efficient for nanotubes p4L EH B3, fRJ. Quantum sieving separates lighter 



molecules from heavier ones by selectively adsorbing heavier molecules. Such 
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selective adsorption can be explained by the higher zero-point energy of the 
light species, which makes their adsorption relatively unfavorable. Quantum 
sieving can be implemented when the adsorbate is effectively confined to a 
one-dimensional (ID) channel or a 2D surface or a small cavity (OD). The 
selectivities observed for adsorption of H 2 isotopes on common substrates 
such as graphite, zeolites, and alumina are low - typically in the range of 1.1 
to 3 |3l], |32], |33|. However, differences in zero-point energies of the adsorbed 
species are expected to be particularly large when molecules are confined in 
very narrow pores. Wang et al. have shown that nanopores with diameters 
wider than 7 A exhibit weak selectivity while smaller nanopores were pre- 
dicted to exhibit large selectivities |24]. Carbon nanotubes typically have 



diameters larger than 7 A, but the interstitial channels are smaller, so that 
they have a pore size and solid-fluid potential which can effectively sieve 
mixtures of H 2 and D 2 . Indeed, Path Integral Monte Carlo calculations at 
variable pressures |25j yielded large isotope selectivity in the IC. Our calcu- 
lations are performed in the limit of low coverage (virtually zero pressure), 
using what we believe to be an improved gas-solid interaction potential. 
The solid-fluid potential is usually modelled as a pair- wise sum of Lennard- 



Jones interactions between adatoms and carbon atoms pR p6| . The potential 



used in our study is the one recently developed by Kostov et al. ||35| . In ad- 
dition to terms present in previous calculations, Rostov's potential includes 
also the interactions of the H 2 quadrupole moment with the local electrostatic 
field of C nanotubes (found from ab initio calculations pi| ) and interactions 
of the H 2 static multipole moments with the image charges induced on the 
surrounding nanotubes. The model yields a strongly confining potential for 
if 2 and D 2 , of depth about 1800 K (~ 160 meV), which alters the rotational 
spectrum and induces a large difference between isotopic zero point energies. 
As we will show, the calculated difference between H 2 and D 2 binding ener- 



gies is consistent with recent isosteric heat data of Wilson et al. ||37|| . This 
large difference, when exponentiated in an Arrhenius expression at low tem- 
perature, yields huge isotope selectivities in the IC and groove, far exceeding 
that on common sorbents (e.g. graphite and zeolites). 

The rotational hindrance of H 2 molecules adsorbed in carbon nanotube 
bundles is a controversial issue. There exist experiments showing no alter- 
ation of the rotational spectrum 3£| and other experiments (with different 
samples) which manifest large changes compared to the free rotation energy 
spectrum [RIJ. The calculations of Kostov et al. indicate that interstitially 



adsorbed H 2 has a significantly hindered rotational motion [35, 36|. The exis- 
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tence of a large rotational barrier leads to large splittings and shifts (relative 
to an orient at ionally invariant potential) of the H 2 and D 2 rotational energy 
levels (see Fig.l). In this work we show that a consequence of this hindered 
rotational motion is a significant spin selectivity in the ICs and grooves. 

Recently, Hathorn et al. showed that quantization of the restricted rota- 
tional motion of H 2 and D 2 confined in SWNTs contributes significantly to 



quantum sieving ||45|| . Their study, however, is substantially different from 
the present work, since they consider H 2 adsorbed within a single nanotube 
and approximate the rotational potential as that of the molecule on the axis 
of the nanotube. In addition, they employed a modified Lennard- Jones po- 
tential with e and a parameters far different from the potential parameters 
assumed in this work. They find significantly hindered rotational spectrum 
for the adsorbed molecule for small nanotubes. Their results are qualitatively 
consistent with our results although we are considering quite different sites. 

The outline of this paper is as follows. In the next section, we calculate the 
ortho-para selectivity as a function of temperature (T), taking into account 
the hindered rotational motion. In Section III, we derive the specific heat 
due to the rotational motion of the adsorbed molecules. Section IV presents 
results for the isotope selectivity of mixtures of H 2 and D 2 . Section V presents 
calculations of the isosteric heats of H 2 and D 2 and compares these with 
the experimental results of Wilson et al. [P7|] . Section VI summarizes and 
discusses this work. 



2 Ortho-para selectivity 

In order to compute the selectivity of component 1 (e.g. ortho) relative to 
component 2 (e.g. para) at low coverage, we use the equilibrium condition 
relating the chemical potential of each component of the 3D vapor (outside 
the nanotube bundles) to that of the ID gas adsorbed in the IC: 

^ ap ° r = tf 3 (2) 

where i stands for each component, ortho and para. We employ statisti- 
cal mechanics to find the chemical potential fi = (dF/dN)T,v, where F = 
—ksTlnQ is the Helmholtz free energy, Q is the canonical partition function 
and ks is Boltzmann's constant. In the low density limit, the partition func- 
tion involves energy contributions from (hindered) rotation, vibration and 
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from translation of the center of mass (CM) in the presence of a transverse 
confining potential. 

In the IC or groove, the molecules are assumed to move freely along the 
axial direction, so the fluid is described by a constant (mean) ID density, p. 
The effect of corrugation on the motion is sensitive to assumptions about 
registry of adjacent tubes [IT) and is neglected here. The CM transverse 



oscillation is described by the one-particle partition function, q cm . The inter- 
nal degrees of freedom consist of the intra-molecular vibrations and hindered 
rotations, described by q mb and q rot , respectively. Therefore: 

Fl c = -Nk B T\n[^-qrqf b (9 l qr t )^M-C l /(k B T))} (3) 

H\ C = Ct- k B T ln(-l-) - k B T In <£" - k B T In qf> - k B T In (4) 

where i stands for para and ortho, A = \J2-kTi 2 / (mk B T), gi is the nuclear spin 
degeneracy and Cj is the interaction energy experienced by molecules due to 
the confining environment. 

The partition function of the coexisting gas outside the nanotube bundle 
has contributions from the 3D translation motion and from internal degrees 
of freedom (vibration and free rotation of the molecules). The free energy 
and the chemical potential are then: 



pvapor = -Nk B T\n[^- q f( gi qtn} (5) 
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— —k B T\a — — - k B T In qf - k B T In ( 9l q{ ree ) (6) 



where n« is the 3D density of ortho or para vapor outside the nanotube 
bundle. 

Experimental data of Williams et al and calculations of Kostov et al 
indicate that the intra-molecular vibrations are essentially the same in the IC 



as in free space |23| , |44fl . The CM oscillations depend only on the molecular 
mass, so they give the same contribution to the ortho and para partition 
functions. Finally, the only relevant contributions to the spin selectivity 
come from the even and odd rotational partition functions, denoted q + and 
g_, respectively: 

Q+= J2 g jm exp(-f3e jm ) (7) 
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where ft = l/(k B T), gj m is the rotational degeneracy (not including the 
nuclear degeneracy), and e jm are the rotational levels. For the odd partition 
function (q~), the summation in Eq.7 is done only over the odd rotational 
levels. In free space, both H 2 and D 2 can be modelled as free quantum rotors 
with gj m = j(j + 1) and e jm = j(j + 1)B. Here j is the rotational quantum 
number and B = h 2 / (21), the rotational constant of H 2 or D 2 . However, the 
different character of the nuclear spin (fermionic for H 2 and bosonic for D 2 ) 
imposes different conditions on the total wavefunction (antisymmetric for 
H 2 and symmetric for D 2 ), so that the nuclear degeneracies corresponding to 
even and odd rotational states are different: g+ 2 = 1, g 1 ! 2 = 3 and g+ 2 = 6, 
g® 2 = 3. By convention, the name ortho is given to the species with larger 
statistical weight, so the even states are called para for H 2 and ortho for D 2 . 
Thus, at high T, the ortho and para molecules exist in the ratio of 3:1 in the 
case of H 2 and 6:3 in the case of D 2 . 

The spin selectivity in the low coverage (pressure) limit involves a ratio of 
densities, which can be found by imposing the equilibrium condition (Eq.2): 

S= . ip - /p ^ IC (8) 

[Jl-- j vapor 

where — and + subscripts correspond to the odd and even rotational states, 
respectively (ortho and para for H 2 , para and ortho for D 2 ). All similar 
factors mentioned previously cancel out in S and, in the end, the low coverage 
selectivity depends on just the even and odd partition functions, q + and g_: 



S= fa' C M°) (9) 
(q va ^ or j q va P or \ 



A qualitative understanding of the spin selectivity can be achieved by 
examining the rotational spectrum. The ortho-para ratio in the IC, in the 
limit of high T, is the same as in the gas phase. But since the energy spectrum 
changes in the IC, the overall partition functions differ, with particularly large 
effects at low T. Fig.l contrasts the rotational spectra of H 2 in the IC (as 
computed by Kostov) and in free space. In the IC, the reference energy is 
the energy the molecule would have if it does not rotate (Ch 2 =-1275 K)p5|. 
Notice that the ground (para) state is shifted down from this value by ~ 200 
K and the first excited state (ortho j = I, rrij = 0) is also shifted down (~ 
60 K), close to the ground state. The D 2 rotational spectrum undergoes a 
similar alteration in the IC, but with a smaller shift of the ground (ortho) 
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state (about 130 K) with respect to a different reference energy ((7d 2 =-1445 
K), and a slightly larger excitation energy to the (J = l,rrij = 0) para state 
(about 85 K)[^]. The shift in the first excited state (ortho for H 2 , para for 
D 2 ) towards the ground state (para for H 2 , ortho for D 2 ) causes the ortho to 
para ratio of H 2 molecules, or para to ortho ratio of D 2 molecules, to increase 
substantially in the IC with respect to that in free space. In the groove site, 
these ratios are still larger than in free space, but smaller, however, than in 
the IC, due to the lesser confinement between the two nanotubes in the groove 
site. At low T, all but the lowest states yield negligible contributions to q rot , 
so that the selectivity depends only on the values of the lowest excitation 
energies in the IC and free space. 

q _ expj-geip - epp)] + 2 exp[-/3(e n - e 00 )] 

3exp[-2/?B] UUj 

where the free space separation is twice the rotational constant B=85 K for 
H 2 . 

The results for the ortho-para selectivity are shown in Fig. 2. We find that 
the spin selectivity is larger than 1 for an extended range of temperatures, 
up to 100 K in the IC and ~ 75 K in grooves. This is consistent with the 
difference between the excitation energies in the IC (groove) and in free space 
which enters Eq.10. As discussed above, the alterations of the rotational 
spectrum favor the H 2 ortho species, whereas for D 2 , the para molecules 
are preferentially adsorbed. Since the H 2 molecules have a larger zero point 
motion than D 2 molecules, they experience more rotational hindrance and 
their spectrum is more altered in the corresponding energy scale. Since the 
reference energy does not count for the spin selectivity, this exhibits larger 
values for H 2 than for D 2 . 



3 Rotational specific heat 

In order to compute the thermodynamic properties of the adsorbate, one 
works with the composite equilibrium partition function: 

g£ = S(2S + l)q + + (S + 1)(25 + l)g_ (11) 

q% = (S + 1)(2S + l)q+ + S(2S + l)g_ (12) 

where S = 1/2 for H 2 and S — 1 for D 2 . We will discuss the H 2 case 
and point out the differences in the case of D 2 . In many circumstances, 
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H 2 is not in thermal equilibrium as regards the relative magnitudes of para 
and ortho components, because the probability of flipping the nuclear spin is 
very small (the lifetime in free space is ~ one year). Therefore, the transition 
probability of a molecule from one nuclear spin state to another is negligible 
during a specific heat experiment. Consequently, the sample may be viewed 
as a non-equilibrium mixture of two independent species, which give additive 
contributions to the heat capacity. 

C n on—eq /pCp foC (13) 

Here f p and f a are the para and ortho molar fractions. In an equilibrium 
mixture these depend on temperature, whereas in a non-equilibrium mixture 
their relative concentrations are determined by the initial conditions. In our 
calculations, we let the para and ortho molar fractions have their high T 
values: f p = 1/4 and f — 3/4. The para and ortho specific heats (in the low 
density limit) are: 

C„ = N^ B ±(T a ±\nq + ,-) (14) 

At low T, only the lowest levels contribute to the specific heat. Let us 
consider only the first two levels: e o and €20 for the para species, and eio and 
e n for the ortho species. Applying the formulae above, the para and ortho 
specific heats are: 

Cp = « >- eM r w "°:';' (i5) 

p (1 -\- e -/?(e20-eoo))2 > 

f ei1 — ein ) e -/ 3 ( e ii- e io) 
C ° = Nk ° (l + 2e-t.,-,..y < 16 > 

Thus, at low T, the para and ortho heats in the IC depend only on the 
excitation energies (e 2 o — eoo) an d (en— eio), respectively. Fig. 3 (a,b) shows a 
comparison between the para, ortho and non-equilibrium specific heats in free 
space and in the IC. Notice that, at low T, in the IC the main contribution 
to the net specific heat comes from C G , whereas in free space it comes from 
C p . This is related to the finite excitation energy (en — eio) between (j=l) 
ortho states in IC, which are degenerate in free space (see fig.l). It can be 
shown that the peak in the heat capacity occurs at T A/(3k B ), where A 
is the first excitation energy present in a specific problem. Therefore, in free 
space the para peak is at about 80 K, while in the IC the ortho peak is at 
about 75 K. 
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In some environments, however, the ortho-para conversion may occur 
quickly (e.g. due to magnetic impurities). In this case, the net rotational 
specific heat is obtained from the composite partition function (Eqs.ll, 12): 

C eq = Nk B -^(T 2 -^\nq rot ) (17) 

At low T, the equilibrium specific heat depends only on the (eio — eoo) exci- 
tation energy: 

C ^ ~ NkB (l + e -^ 10 -coo))2 ^ 

The equilibrium specific heats of H 2 and D 2 in the IC are shown in Fig.4. 
In the case of H 2 , there is a substantial ortho-para conversion peak at about 
20 K, followed by a gentle rise above 75 K, corresponding to ortho-ortho 
excitation. However, in the case of D 2 , the ortho-para conversion appears 
only as a small bump at about 25 K, since the molar weights of para and 
ortho species are now different. At high T, the rotational specific heat per 
molecule goes to k B) the classical result for a free rotor. 



4 Isotope slectivity 

Let us now consider a mixture of H 2 and D 2 . The equilibrium condition 
between the ID gas adsorbed in the ICs or grooves and the coexisting vapor 
outside the nanotube bundles is given again by Eq.2: ii\ c = fiJ apor , where 
now % stands for H 2 or D 2 . The partition functions and the chemical po- 
tentials of each species consist of the same factors, but their meanings are 
different. For example, the rotational partition function is given now by the 
composite partition function (_Egs.ll, 12) of each isotope. Since the H 2 and 
D 2 species have different masses, their zero point energies will be different, 
and these will give a nontrivial contribution to the selectivity (as will dif- 
ferences in the rotational energies). The transverse CM oscillations of the 
molecule in the IC or grooveare treated as excitations of a 2D harmonic 
oscillator. The corresponding partition function is: 

te = Ee-^» + V,, ]2 = _^__ (19) 

where oo is the frequency of CM oscillations, found from the force constant. 
The rotational hindrance is expressed in terms of the rotational energy shifts, 
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e m j and the additive constants (C) in the expansion of the rotational part of 
the interaction potential |35 . This additive (reference) energy is also diferent 



for D2 and if 2- The rotational hindrance leads to a shift of the binding energy 
of each species: Ef, = — {C + Tiuj + eoo)- 

The equilibrium condition for each species (i = D21H2) yields the ratio 
of IC (p^ and vapor (n^) densities for each species, which is then used to find 
the isotope selectivity. 



Pi _ , 2 exp(-Phuj) exp(-pC)q^ t 
— ~ \~ / ot. . w>. JFTe K Z[J ) 



Hi 1 (1 - exp (-(3huJi)) 
The isotope selectivity is: 

S = (21) 

n D .Jn H2 

Eq.20 implies that at low T the main contribution to selectivity comes from 
the difference between the binding energies 

S « i exp \-(3{C D2 + hu D2 + e D 2 - C H2 - hco H2 - eg*)} = \exp[(3(E^ - E^)} 

(22) 

Fig. 5 displays the results of our calculations for isotopic selectivity in the 
ICs and grooves of an (18,0) nanotube bundle, and on graphite (for compar- 
ison). All three surfaces favor D2. However, the isotope selectivity in the IC 
and groove is much larger than that on graphite and within the tubes . 
Even at 300 K the selectivity in the IC and groove is substantial, due to the 
large difference between the D2 and H2 binding energies. This large differ- 
ence comes primarily from the zero point energy (but the difference between 
the rotational energies also contributes). Table 1 presents our results for the 
different contributions to the the binding energies. The difference between 
D 2 and H 2 binding energies is ~ 400 K in the IC and ~ 170 K in the groove, 
much larger than the difference of ~ 35 K on graphite |BS| . When exponenti- 



ated, these differences lead to the very large selectivities found in the IC and 
groove. We may compare our results with previous calculations. Challa et 
al. have studied isotope selectivity at variable pressure in the interstices of 
(10,10) nanotube bundle, using Ideal Adsorbed Solution Theory(IAST) and 
Path Integral Grand Canonical Monte Carlo simulations (PI-GCMC) ||25|| . At 
low pressures, the IAST and PI-GCMC results converge to the analytic zero 
pressure results. In their calculation, only the transverse zero-point energy 
is considered, i.e. rotational motion is ignored. The selectivity calculated 
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by them at 20 K and very low pressure (10 -13 atm) is about 700. (10,10) 
and (18,0) nanotubes have similar diameters (13.6 and 13.8, respectively), so 
the interstices of a closed packed nanotube bundle are similar, too (~ 6 A 
diameter). Recently, Hathorn et al. [[45] have estimated the quantum effects 
due to the restriction of the rotational motion of H 2 and D 2 adsorbed within 
a single nanotube. Their results reveal significantly enhanced rotational sep- 
aration factor (selectivity), which at 20 K is of order 10 3 . The zero pressure 
isotope selectivity calculated by us at 20 K, in the interstices of bundles of 
(18,0) carbon nanotubes is of order 10 9 . Our results differ to such a large 
degree since we considered both zero point and rotational contributions to 
selectivity and we used a slightly deeper, solid-fluid potential. The stronger 
the confinement, the larger the difference between the isotope zero point en- 
ergies and between the rotational energies. Previous calculations inside (5,5) 
nanotubes at 20 K yielded much lower isotope selectivities PBJ. At low pres- 



sure, the computed selectivity of tritium to H 2 is 23, and that of tritium to 
D 2 is 1.7; hence Sq(D 2 / H 2 ) = 13. Thus, even if the inner channels of (5,5) 
nanotubes have about the same radius (3.3 A) as the ICs of (18,0) nanotube 
bundles (3 A), the computed selectivity in the IC is much larger due to the 
different arrangement of the carbon atoms in the IC, which yields greater 
confinement. 



5 Isosteric heat 

The isosteric heat is defined as: 

where P is the pressure of the vapor outside the nanotube bundles. In our 
model, at low density, the dependence of P on T can be found from the 
equilibrium condition (Eq.2) using the ideal gas formula: n = f3P. Finally, 
the isosteric heat is: 

Q st = 2k B T -Hu-C + i-) [In ^ZJy ] ^ 

At low T, the excitations can be neglected and the term involving the 
derivative with respect to (3 goes to — e o — (eio — e o)e~ /3 ^ 10_<:00 ^ + 2Be~ 2/3B , 
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where eoo and eio are the first two rotational energy levels. Therefore, at low 
T, the isosteric heat is a measure of the binding energy: 



Q st = E b + 2k B T (25) 

since Ef, = — (C + Tiu + eoo)- 

The isosteric heat can be found experimentally by taking the difference 
between two nearby isotherms. Recently, there have been reported experi- 
ments of adsorption of D 2 and H 2 in carbon nanotube bundles [33], which 
show isosteric heats in nanotube bundles to be a factor of 1.5 (for H2) to 
1.8 (for D2) larger than those on graphite. This means that the D2 and H2 
binding energies are almost twice as large as those on graphite. Moreover, 
the difference between isotope isosteric heats, which again is related to the 
difference between isotope binding energies, was found to be about 200 K, 
much greater than the difference of about 35 K on graphite f39 ]. This is 



a mark of the greater confinement of molecules adsorbed in nanotube bun- 
dles than on the surface of graphite; the confinement leads to the enhanced 
separability of isotopes found in the previous section. 

We performed calculations of the isotope isosteric heats in the low den- 
sity limit (Eq.24). Table 2 and Figs. 6(a),(b) compare the results of the 
calculated isosteric heats and the low coverage experimental data at T=85 
K. The experimental procedure of Wilson et al does not reveal where in 
the nanotube bundle the adsorption occurs. Since the experimental isosteric 
heats are much larger than on graphite, one assumes that the adsorption 
environment should be either the IC or the groove site. Our calculations of 
the isosteric heats in both IC and groove are lower than the experimental 
results. The calculated difference between D 2 and H 2 isosteric heats (374 
K) in the IC is much larger than the one found experimentally, but the dif- 
ference (137 K) in the groove is closer to the experimental difference (200 
K). In addition, the calculated binding energies(see Table I) for D 2 and H 2 
in the grooves are consistent with the experimental results and show that 
the adsorption on the grooves is much more likely than on the IC at low T. 
Assuming the accuracy of our approximations, we believe that our results 
are compatible with an adsorption in the groove, rather than the IC. This is 
also experimentally plausible, since the external grooves are more accessible 
for adsorption than the ICs. Indeed, the ICs may be completely blocked, 
preventing any adsorption there. 
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6 Summary and conclusions 



We have investigated the spin (ortho versus para) and the isotope (D 2 versus 
H 2 ) selectivity at low coverage in the IC and groove channel of (18,0) nan- 
otube bundles. The rotational hindrance of molecules in the IC and groove 
site induces shifts and splittings in the rotational spectrum, enabling the 
nanotube bundles to preferentially adsorb H 2 ortho molecules and D 2 para 
molecules. Our calculations show substantial spin selectivity for tempera- 
tures up to ~ 100 K in the IC and up to ~ 75 K in the groove. At low 
temperatures, the alterations of the rotational spectrum in the IC induces 
new features in the ortho heat capacity, as the ortho-ortho peak. The H 2 
equilibrium heat capacity exhibits a distinctive ortho-para conversion peak 
at ~ 20 K. The different features of non-equilibrium and equilibrium heat 
capacities may be a way to check the existence of the ortho-para conversion 
of H 2 and D 2 adsorbed in nanotube bundles. In the case of the isotope selec- 



tivity in nanopores, zero-point motion favors the heavier isotope p5f, and the 



rotational motion enhances the preferential adsorption of the heavier isotope 
p5| , [R| . We find substantial isotope selectivity even at temperatures as high 
as 300 K in the IC and groove. At 20 K, our calculation at low coverage 
(consequently zero pressure) yield isotope selectivities of order 10 9 in the IC, 
orders of magnitude larger than the low pressure limit found in the calcula- 
tions of Wang et al. and Challa et ai. Our much larger isotope selectivity 
is a consequence of the stronger confinement in the IC, which yields larger 
zero point and rotational energies. The difference between the D 2 and H 2 
isosteric heats in the groove (~ 135 K) is close to the value found experi- 
mentally (~ 200 K). We conclude that carbon nanotube bundles are ideal 
surfaces for spin and isotope separation, providing a potential technology of 
quantum sieving. 

We acknowledge stimulating discussions with J. K. Johnson, M. M. Calbi, 
T. Wilson, M. J. Bojan, O. Vilches, P. Sokol, D. Narehood and D. Stojkovic. 
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Table 1: The zero point energies (hu), rotational energy references (C) and 
shifts (e 00 ) and binding energies (E b ) (expressed in Kelvin). D 2 — H 2 means 
the difference between the corresponding values. 





hw IC 


^groove 




^groove 


IC 

fc oo 


groove 

e oo 




■^groove 


j^graphit 


D 2 


762 


266 


-1445 


-970 


-133 


-62 


816 


766 


517 


H 2 


1077 


370 


-1278 


-880 


-206 


-93 


407 


603 


482 


D 2 -H 2 


-315 


-114 


-166 


-90 


73 


31 


409 


163 


35 



Table 2: The isotope isosteric heats (Q) in the IC, groove and experimental 
results (expressed in Kelvin) 





QIC 


^.groove 


Qexp 


D 2 


991 


901 


1100±90 


H 2 


617 


764 


900±70 


D 2 -H 2 


374 


137 


200±110 
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Figure captions 

Fig. 1. H 2 rotational spectrum in free space and in IC. The number of 
lines in a level represents the degeneracy of that energy level, o and p mean 
ortho and para, respectively. 

Fig. 2. D 2 and H 2 spin selectivity in IC and groove channel as a function 
of T 

Fig. 3. Ortho (dashed line), para (dotted line) and non-equilibrium (full 
line) specific heat (without spin equilibration) of H 2 molecules (a) in IC, (b) 
in free space 

Fig. 4. The equilibrium rotational specific heat of H 2 (full line) and of 
D 2 (dashed line) 

Fig. 5. Isotope selectivity in the IC (full line), groove (dashed line) and 
on graphite (dotted line). 

Fig. 6. The calculated H 2 (dotted line), D 2 (dashed line) and the differ- 
ence D 2 — H 2 (full line) isosteric heats (a) in the IC and (b) in the groove. 
The experimental values at 85 K are shown in symbols: circle for D 2 , square 
for H 2 and triangle for the difference D 2 — H 2 . 
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